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Spontaneous CP-violation in the strong interaction is analyzed at 9 = n within the framework 
of the two-flavor NJL model. It is found that the occurrence of spontaneous CP-violation at 9 = 
7r depends on the strength of the 't Hooft determinant interaction, which describes the effect of 
instanton interactions. The dependence of the phase structure, and in particular of the CP-violating 
phase, on the quark masses, temperature, baryon and isospin chemical potential is examined in 
detail. When available a comparison to earlier results from chiral perturbation theory is made. From 
our results we conclude that spontaneous CP-violation in the strong interaction is an inherently low- 
energy phenomenon. In all cases we find agreement with the Vafa-Witten theorem, also at nonzero 
density and temperature. Meson masses and mixing in the CP-violating phase display some unusual 
features as a function of instanton interaction strength. A modification of the condition for charged 
pion condensation at nonzero isospin chemical potential and a novel phase of charged ao mesons are 
discussed. 

QJ • PACS numbers: 12.39.-x,11.30.Er,11.30.Rd 
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I. INTRODUCTION 



r— i | The possibility of CP-violation in the strong interaction has been studied extensively. The QCD Lagrangian 

i' naturally incorporates a #-term Cg = 32 ^. 2 FF, which can lead to CP-violation due to instanton contributions. 

in: 



32tt 2 

Only for 9 = (mod tt) the Lagrangian is CP-conserving. In nature the value of 9 is extremely close to zero, as 
has been concluded from pseudoscalar mass ratios and the neutron electric dipole moment [H, 0, Q . This suggests 
that 9 = 0, but the lack of a satisfactory explanation of why this should be the case is commonly referred to as 
the strong CP problem. 

At 9 = no explicit CP-violation is present in the Lagrangian, but in addition, the well-known Vafa-Wittcn 
theorem states that spontaneous parity violation in QCD at 9 = docs not arise Q ■ This rules out the possibility 
Q\ ' that (FF) ^ at 9 = and implies that the QCD ground state for nonzero 9 must have higher energy than at 
, 9 = 0. The situation is different at 9 = tt, when the Lagrangian is also explicitly CP-conserving. In this case 
^>0 spontaneous CP-violation could arise, as was first pointed out by Dashen [f|. There are then two degenerate 
CP-violating vacua, which differ by a CP transformation from each other. 

This possibility of spontaneous CP-violation is one of the reasons why people have studied the (9-dcpcndence of 
the strong interactions. However, to study this in full QCD is very difficult due to the nonperturbative nature of the 
0-tcrm. Even in lattice QCD, studies are limited to small 9, because of the problem of how to deal with complex 
phases. Therefore, the ^-dependence of the strong interaction and Dashen's phenomenon have been studied 
extensively using low energy effective theories, such as chiral perturbation theory [1, 0, @, H, GU EH, E2, EH , or by 



using specific models, such as the NJL model |14[. In a quark model, like the NJL model, the effects of instantons 



and the #-term are incorporated via an effective interaction, the 't Hooft determinant interaction [1J5|, |16fl . In chiral 
perturbation theory these effects can be included in a similar way via a log determinant interaction^, 0] . Here 
we shall denote the strength of the latter interaction by a. Whether or not the theory exhibits spontaneous CP- 
violation at 9 = 7r depends on a and on the values of the quark masses. Two limiting cases were discussed in the 
literature. Ref. Q, which considers lowest-order chiral perturbation theory, states that when a/N is nonzero but 
much smaller than the quark masses, the theory always exhibits spontaneous CP-violation at 9 = tt, independent 
of the values of the quark masses and number of flavors N. The opposite case 0, i.e. when the masses 
of the quarks are much smaller than a/N, leads to different results. In this case it does depend on the values 
of the quark masses. In the two-flavor case for a/N —> oo (which means no t] meson is included), spontaneous 
CP-violation only occurs for degenerate quark masses (in that case actually for all finite values of a/N). For 
finite a/N 3> m u , spontaneous CP-violation also occurs for nondegenerate quark masses in a finite interval of 
nid/m u around 1, as was shown in Ref. Q. In the three-flavor case, a region exists in the (m u , TOd)-plane where 
the theory spontaneously violates CP invariance [ll[ , as shown in Fig. [TJ The asymptotes depend on the value of 
the strange quark mass. 

In first-order chiral perturbation theory there are only a few parameters, namely the quark masses, the pion 
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FIG. 1: The (m«,m<j) phase diagram of three-flavor chiral perturbation theory, with 6 = tt and a/N — > oo, cf. Ref. [ill ). 

decay constant and the strength of the determinant interaction. It is therefore interesting to study CP-violation 
in a somewhat richer situation, such as chiral perturbation theory beyond leading order, which has been studied in 
Ref. Q for a/N — > oo. Here we will make a comprehensive study of the ^-dependence and especially spontaneous 
CP-violation at 9 = tt within the framework of the two-flavor NJL model in the mean-field approximation. We 
will study the dependence on the effective instanton interaction strength (denoted by c in this case), not only in 
the two limiting cases, but for all possible values. We find that there is a critical value of the interaction strength 
above which spontaneous CP violations occurs and which depends linearly on the quark masses, as expected from 
axial anomaly considerations. 

As will be discussed, the two-flavor NJL model allows for spontaneous CP-violation also for nondegenerate quark 
masses. We find a region in the (m u , m^-plane very similar to the three-flavor lowest-order chiral perturbation 
theory result shown in Fig. [TJ However, for the two-flavor NJL model the asymptotes are determined by the 
strength of the instanton induced interaction. 

We will also study the influence of nonzero temperature and baryon and isospin chemical potential. It has 
been suggested that in those cases the Vafa-Witten theorem may no longer apply (see for instance Ref. [l7| for 
some explicit arguments, but also Ref. [H[ for counterarguments). But even if it does apply, spontaneous CP- 
violation at finite temperature or baryon chemical potential (through meta-stable states) has been considered in 
the literature [ljj, H3, Hi], H2] and possible experimental signatures in heavy ion collisions have been put forward 
[23l . l24l . I25I . This would also be relevant in the early universe, when possibly 9 was nonzero and later relaxed 
to zero, for example via a Peccci-Quinn-like mechanism [27], [H, [29[. We therefore wish to check the Vafa-Witten 
theorem and the possible presence of CP-violating local minima in the NJL model at finite temperature and 
density. 

Spontaneous CP-violation at 9 = tt within the two-flavor NJL model including temperature dependence has 
been considered before in Ref. [lj], but only for a very limited range of quark masses: \m u ± < 6 MeV and 
without chemical potentials. 

In Ref. [IH the phase diagram as a function of 9 and isospin chemical potential has been investigated within 
first-order chiral perturbation theory for two flavors and effectively a/N — > 00 (due to the absence of the 77 meson). 
We will compare this to our results at nonzero isospin chemical potential, where a modification of the pattern of 
charged pion condensation is observed at 9 = tt. 

This article is organized as follows. First the NJL model is briefly introduced to set the notation, then we discuss 
the effect of chiral transformations on the theory, which is relevant for the calculation of the effective potential 
and for a comparison to earlier results from the literature. We continue with a discussion of the ^-dependence 
of the ground-state, including temperature effects and nonzero baryon and isospin chemical potential. Also we 
discuss the c-dependence of the meson masses and mixing in the CP-violating phase. We end with conclusions 
and a further discussion of the results. 



II. THE NJL MODEL 

The Nambu-Jona-Lasinio (NJL) model, introduced in Refs. [13, [111, is a low energy effective theory that contains 
four-point interactions between the quarks. In this article the following form of the NJL-model is used 

C = ip (i^d/j, + ip - C M + C qq + £ de t, (1) 
where the mass term of the Lagrangian is 

Cm = ^MoV>, (2) 

and [i = (fj. Ul [id) denotes the chemical potential. Furthermore, 

C qq = G x [$\ a ipf + {^Xail^f] , (3) 

is the attractive part of the qq channel of the Fierz transformed color current-current interaction [32| and 

£det - G 2 e je dct(^ i )+h.c, (4) 

is the 't Hooft determinant interaction which depends on the QCD vacuum angle 9. Often G\ and Gi are taken 
equal, which at 9 = means the low energy spectrum consists of a and tz fields only. We will restrict to the two 
flavor case, using A a with a = 0, 3 as generators of U(2). We will not take into account diquark interactions, 
and therefore do not consider color superconductivity that is expected to arise at high baryon chemical potential 
and low temperatures. We choose an appropriate basis of quark fields, such that the mass-matrix Mo is diagonal, 
i.e., 

(*?JL)- < 5 > 

The symmetry structure of the NJL model is very similar to that of QCD. In the absence of quark masses 
and the instanton interaction there is a global SU(3) c xU(2)ftXU(2)L-symmetry. The instanton interaction 
breaks it to SU(3) c xSU(2)i,xSU(2).rxU(1)b- For nonzero, but equal quark masses this symmetry is reduced to 
SU(3) c xSU(2)y xU(1)b. For unequal quark masses and chemical potentials one is left with SU(3) c xU(1)b xU(1)j, 
where B and / stand respectively for baryon number and isospin. 



III. CHIRAL TRANSFORMATIONS AND NEGATIVE QUARK MASS 

It is well known that a theory with 9 = tt can be related to a theory with a negative quark mass. Since this 
sometimes leads to confusion concerning the terminology used for the meson spectrum, we will elaborate on this 
relation in this section. 



A. QCD 

We start with the QCD partition function including the 9 term: 



where 



,2 



Cg = -J^WF^. (7) 

OZ7T 

The fermion measure of gauge theories is not invariant under chiral transformations [33j . which can be used to 
remove Cg. Since the mass term is not invariant under chiral transformations, a 9 dependence then appears in 
the mass term. One obtains 

Z = J ViP'VftVA e'-f^'^, (8) 
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where the (9-dcpcndence resides in the mass term. Although the physical results one obtains using the transformed 
expression will be equivalent, one has to be careful when evaluating vacuum expectation values. We define vacuum 
expectation values of an operator O = 0(V>, V>) in terms of the original fields (the one of the Lagrangian ([7])) and 
in terms of the transformed or "primed" fields as follows: 

{0) g = fV^V^VAOiip,^) e 1 ^ 4 ^ 00 ^, 

(0') g = fVip'-DfiVA O' e i f d4xL '<* CD . (9) 

Clearly, the condensates {0) g and (0') g differ for 8^0 and are related by a ^-dependent transformation. For 
instance, 

(^)e = I T>i'Vi)VA -ipip e'/d^QCD+c*] ^ f jy^'V^'VA flip' j S dixC ' Qcn = (^'ip') . (10) 



When discussing a vacuum expectation value like (4 , ' t P}g it has to be accompanied by a statement about which 
Lagrangian one is using. 

In what follows we will select the chiral transformation that only affects the up quark: 

u L = e~ ie/2 u' Ll 

u R = e i9 / 2 u' R . (11) 
This removes Cg from the Lagrangian and the up-quark mass term changes according to 

um u u = u [m u cos 9 + m u i^ sin 9] u' . (12) 
For 9 = tt a negative up-quark mass results. In addition, 

(uu) g — (u'u') g cos 9 + (u'ij5u') g sm8. (13) 
Below we will use the following notation for the meson condensates: 

(a) = $\ o1 b) , (ao) = (W) 

(rj) = (^75^) , (w) = (tMhsV') ■ (14) 

These condensates transform according to: 

(a) = \ (cos 9+1) (a'} + \ (cos 6 - 1) <a° ') + \ sinfl (r/> + \ sin0 (tt ') , 

(a Q ) = cos - (a } + sm - (tt* ) , 

(a°) = \ (cos 6-1) + ! (cos 0+1) (a° ') + \ sin 6 (rf) + \ sin 6 (tt ') , 

(V) = I (cos 9+1) ( ? /) + i(cos0-l) (O - isinfl (</) - isinfl (ag 7 ) , 



(tt*) = cos^ (Tr^-sin^ (4') 



(tt ) = l( C0S e-l) (r,') + ^( C o S e + l) (7r°')-l S me (a')-l S me (a° '). (15) 

Therefore, one has to be careful assigning the names ir° and 77 to the condensates after doing a chiral transforma- 
tion. For example, Ref. [lij discusses a (7r°)-condensate using a Lagrangian without 9 term, but with negative up 
or down quark mass. This corresponds to an (77) -condensate using a Lagrangian with positive quark masses and 
a (9-term with 9 = tt. We emphasize that these transformations are just a matter of consistently naming mesons 
and vev's, but this is nevertheless important for the comparison of quantities from different calculations. 
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B. NJL-model 



As the NJL-model is not a gauge-theory, the fermion measure is invariant under chiral transformations. But 
now the Lagrangian contains two terms that are not invariant under chiral transformations, the mass-term and the 
determinant interaction. The latter is 0-depcndent. Like for QCD, this (9-dcpendence can be put in the up-quark 
mass term using a chiral transformation. So the analysis for the NJL-model is similar to the analysis for QCD, 
but instead of a noninvariant measure we have a noninvariant effective interaction. 

The calculation of the ground-state of the NJL-model is more conveniently done with the ^-dependence in the 
up-quark mass term, i.e. we use 

C'm = u' R m u e~ l6 u' L + d' R mdd' L + h.c, 

£'det - G 2 det (V^Vl) + h.c. (16) 

Therefore, below we will calculate the effective potential using the transformed (primed) fields, but discuss the 
ground state phase structure solely in terms of the condensates in terms of the original fields. Only in the latter 
case the SU(2)y symmetry among the three pions (and among the ao-mesons) is manifest when we consider 
degenerate quark masses for instance. 

Because we want to investigate the effects of instantons on the vacuum, we are interested in the dependence 
on the strength of the determinant interaction, which is the effective instanton interaction. Frank et al. [34j have 
investigated the effects of this interaction at 8 = 0, in particular flavor-mixing effects, on the QCD phase diagram, 
by choosing the following expressions for G\ and G 2 (where our c is their a) 

Gi = (1 - c)G , G 2 = cG . (17) 

In this way, the strength of the instanton interaction is controlled by the parameter c, while the value for the quark 
condensate at 6 — (which is determined by the combination G\ + G2) is kept fixed. As mentioned, for G\ = G 2 , 
or cquivalcntly c = i, only the a and 7T mesons are present. For our numerical studies we will use the following 
values for the parameters unless stated otherwise: m u = m& = 6 MeV in case of degenerate quark masses, a 
three-dimensional momentum cut-off A = 590 McV/c and Go A 2 = 2.435. This corresponds [H[ to a pion mass 
of 140.2 MeV, a pion decay constant of 92.6 MeV and finally, a quark condensate (uu) = (dd) = (—241.5 MeV) 3 . 
These values are in reasonable agreement with experimental determinations. 



IV. CALCULATION OF THE GROUND-STATE 



To obtain the ground-state of the theory, we introduce 8 real condensates as follows 

a' = -2(G 1 + G 2 )(a'), 

a' = -2(G 1 -G 2 )(a' ), 

& = -2(Gi - G 2 ) ( V ') , 

= -2(G 1 + G 2 )(n'). (18) 

All quantities in this section refer to the primed fields, but for notational convenience we will drop the primes from 
now on in this section only. Results presented in the subsequent sections will refer exclusively to the unprimed 
quantities. 

It will be assumed that all condensates are space-time independent. A Hubbard-Stratonovich transformation 
eliminates the four-point quark interactions, such that the Lagrangian becomes quadratic in the quark fields 
and the integration over these fields is straightforward to perform. One obtains the following expression for the 
thermal effective potential V in the mean-field approximation [35| 

v = W^h) + W^G 2 -)~ T J W logdctK (19) 

po = (2n+l)7ri 

where if is a matrix in flavor and Dirac space, 



K = l f <g> (i^oPo + HPi) ~ M ® 7o - M (20) 
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is the inverse quark propagator, and 

M = m u (cos6 X u <g> l d + sin^ X u <E> i^ 5 ) + m d X d <E> l d + a a X a ® l d + faX a ® «75, (21) 

with A„ = (Ai + A 3 )/2 and X d = (Ai - A 3 )/2. 

The values of the condensates are found by minimizing the effective potential with respect to these condensates. 
By exploiting U(l) flavor symmetry one only has to study the condensates ao, oti, a$, fa, fa, an d fa- In Ref. [IH 
the fa and fa condensates have been ignored based on the Vafa-Witten theorem. As wc wish to check the validity 
of this theorem at finite temperature and density in our model calculation, we do take these condensates into 
account. 

In order to calculate the effective potential, it is convenient to multiply K with If <E> 70 which leaves the 
determinant invariant and yields a new matrix K with zpo's on the diagonal. The determinant of K can be 
calculated as dct K = Y\ 8 i=1 {Xi — ipo), where Aj are the eigenvalues of K with p = 0. After performing the sum 
over the Matsubara frequencies, we obtain 

r M , ^ n „„A, , .-l^lMl (22) 



T logdetA' = ^ 

p = {2n+l)-nT i=l 



^+Tlog(l + e -' A *'/ T 



Finally we need to integrate over the three-momenta p up to the ultraviolet cutoff A to determine the effective 
potential. 



Minimizing V implies solving the equations 



dV 

aST ' < 23 > 



where x = {ao,a\,a^,fa,fa,fa}. The derivatives of the effective potential can be calculated from the expres- 
sion [3^ 

E logdetA = ^6 tJ (l- p|Ai|/ 2 )sgn(A a ), (24) 

3 p = (2n+l)7rT i=l V / 

where hj = (WdK(p = 0)/dxjU) . Here U is a unitary matrix which contains in the i-th column the normalized 

V / ii 

eigenvector of K with eigenvalue A^. Again, one has to integrate over p to obtain the complete derivative. Since 
this calculation does not use the finite distance method, the derivatives can be determined very accurately. Also, 
it is very efficient as one needs the eigenvalues of K anyway in order to calculate the effective potential. 

When a solution to Eq. (|23[) has been found, it has to be checked whether the solution is indeed a minimum and 
not a maximum or saddle-point. This is checked by verifying that the Hessian of the solution only has positive 
eigenvalues. If more than one minimum is found, the one with the lowest value is chosen. Also the continuity of 
the effective potential is checked. 

The speed of the calculation mainly depends on how fast the eigenvalues of K can be calculated. To speed up 
the evaluation of the calculation of the eigenvalues, one can make use of the fact that the determinant of K is 
invariant under the interchanging of rows and columns. This can be used to bring A to a block-diagonal form of 
two 4 x 4-matrices. This reduces the computing time to determine the eigenvalues with a factor of four as the 
time to numerically calculate the eigenvalues scales cubically with the dimension of the matrix. 

Another way of improving the speed of the calculation is to choose p to lie along the z-direction, exploiting the 
fact that det K does not depend on the direction of p. 

One final remark we have to make regarding Eq. is the fact that in order for the effective potential to have a 
minimum at finite values for the condensates, the coupling G2 has to satisfy —G\ < G% < G\, and correspondingly, 
— \ < c < |. From Eq. (dJ we can see that a negative value for G2 corresponds to shifting 9 — > 9 + n, implying 
that the minimum of the theory will be at 9 = tt, in violation of the Vafa-Witten theorem at zero temperature 
and density. Therefore, we will restrict to < c < |. The case c = | is special, because then only the a' and it' 
fields are present in the theory, which means at 9 = the a and it mesons and at 9 = ir the ?y and ao mesons. 



V. THE GROUND-STATE OF THE NJL MODEL 



This section deals with our results for the ground-state of the NJL model. First we discuss the ^-dependence 
of the condensates and the effective potential. It turns out that for nonzero c, two different situations can be 
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distinguished: below a certain critical c value (c cr it) no spontaneous CP violation takes place at 9 = tt, whereas for 
c larger than this critical value it docs take place. The value of this c cr it depends on the values of the quark masses. 
In Fig. [5] we show the phase diagram at 9 = tt in the (c, ?n)-plane for degenerate quark masses m u = = m, 
two phases can be distinguished, 

1. (cr) 7^ 0, the ordinary chiral condensate. 

2. (a) ^ 0, (?/) ^ 0, the CP-violating phase. 

The phase transition corresponds to c cl -it and is of second order. A linear relation exists between the quark mass 
and Ccrit- 
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FIG. 2: The (c, m) phase diagram at 8 — it. 



A. The ^-dependence of the vacuum 



When the determinant interaction is turned off, there is no ^-dependence. In terms of the unprimcd fields, only 
the (cr) condensate is nonzero. 

In Fig. [3] we show the ^-dependence of the various condensates for the case c = 0.005, which for our choice 
m u = irid = 6MeV is below c cr it « 0.008. As can be seen, no spontaneous CP violation occurs, since (77} = 
at 9 = tt. Explicit CP violation for other values of 9 does occur, as expected. In this figure the condensates are 
normalized with respect to (cr) at 9 = 0. Both (n) and (do) are zero for all 9 and this remains true for c above 
c cr it for degenerate quark masses. 

Fig. H] shows the case of c = 0.2. Spontaneous CP violation is clearly visible, as (rj) is nonzero at 6 = tt. As 
can be seen two degenerate vacua then exist, with opposite signs for (rj). These two degenerate vacua differ by 
a CP transformation. This is known as Dashen's phenomenon Q and is also apparent from the 0-dependence 
of the effective potential. In Fig. \E\ we show the effective potential as a function of 9 normalized to its value 
at 9 = 0, for the two cases c = 0.005 and c = 0.2. In both cases, the minimum of the effective potential is at 
9 = 0, in agreement with the Vafa-Witten theorem. Furthermore, it can be seen that the case with spontaneous 
CP-violation has a cusp at 9 — tt, and therefore a left and a right derivative which differ by a sign. Due to the 
axial anomaly, the 0-dcrivative of the effective potential is proportional to (rj) . This explains the occurrence of 
two values for the rj condensate. 



B. Phase structure at 6 = tt 



In this section we concentrate further on the case 9 = tt. We will start with a discussion of the mass-dependence 
of the ground-state. From Ref. pd| we know that in three-flavor chiral perturbation theory a region exists in the 
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FIG. 4: The ^-dependence of the normalized condensates, with c = 0.2 > Ccrit- 

(m u , md)-plane where CP is spontaneously violated, cf. Fig.[TJ In that case the shape of the CP violating region 
depends on the strange quark mass. In the present case it depends on the choice of c. 

In Fig. E] we show the phase diagram of the NJL model at 6 = n with c = 0.4 in the (m u , m<j) -plane. Four- 
phases can be distinguished 

1. (a) < 0, (og> < 

2. (a) < 0, (eg) > 

3. (a) < 0, (og> < 0, (r,) ? 0, (tt ) + 

4. (a) < 0, (a°) > 0, (r?) ^ 0, (tt ) ^ 

In phases [3] and [3] two degenerate vacua exist with opposite signs for both (77) and (tt°). The phase transitions 



9 




^ 1 ' 1 1 1 1 1 ' 1 

50 100 150 200 250 300 350 400 450 

m u (MeV) 

FIG. 6: The (m u ,md) phase diagram of the NJL at 8 = n with c = 0.4. The dashed lines denote second order phase 
transitions and the dotted line a crossover. 



between the CP-conserving phases \T\ and [5] to the CP-violating phases [3] and 0] are second order. The phases \T\ 
and [2] only differ in the sign for the {a[])-condcnsate, the same holds for the phases |3j and |4] The phase transition 
between the phases [3] and 2] is a crossover, as is the case for the phase transition between phase [T] and [5] for large 
m u and m^. Exactly at the crossover, (aj]} vanishes and in the CP-violating region the same applies to (7r°), 
but not to (77) . The fact that ag-condensation (and 7r°-condensation in the CP-violating region) occurs when the 
masses are not equal simply reflects the explicit breaking of SU(2)y which occurs for nondegenerate quark masses. 

The shape of the CP-violating region is determined by the asymptotes, which arc proportional to c. We conclude 
that in contrast to two-flavor lowest-order chiral perturbation theory (the case of m s — > 00 in Ref . [ll[ , such that 
the asymptotes are moved to m u — rrid = 00), the NJL model does have a spontaneous CP-violating phase for 
two nondegenerate quark flavors. This is in accordance with the chiral perturbation theory analysis of Ref. 0] in 
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the large N c limit and for finite a/N ~3> m u ,md- 



VI. FINITE TEMPERATURE AND BARYON CHEMICAL POTENTIAL 



In this section wc turn to the changes in the phase structure at nonzero temperature and density. Ref. fl4j | 
states that the CP-violating phase at = ir does not exist at high temperatures, i.e. a critical temperature exists 
above which the CP-violating condensates are zero. Ref. [l4| only considered the case c = 0.5. Here we generalize 
their results to other c values. In Fig. [7] the (T, c) phase diagram is shown for degenerate quark masses. The 
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FIG. 7: The (T, c) phase diagram of the NJL model at 
6 — tv. The dotted line represents a crossover, which is 
defined by the inflection points d 2 (a) /dT 2 — 0. 



FIG. 8: The (fj,B,c) phase diagram of the NJL model 
at 9 = 7r. 



following three phases arise 

1. (a) 7^ 0, the ordinary chiral condensate. 

2. (a) ^ 0, (77) ^ 0, the CP-violating phase. 

3. (a) p» 0, the (almost) chiral symmetry restored phase. 

The phase structure at T = can be understood from Fig. for degenerate quark masses the two phases are 
encountered on its diagonal. The phase transition occurs at that particular value of m u = for which c = 0.4 
is the critical c. The phase transition between phases 1 and 2 is of second order for all temperatures. For 
nondegenerate quark masses these two phases would be phases [T] and [3] or [5] and H] of Fig. \§\ depending on whether 
m u is larger or smaller than m,d, respectively. In that case two second order phase transitions are present. 

Above a certain temperature one observes in Fig. [7] an approximate restoration of chiral symmetry (phase 3). 
Note that the chiral symmetry is not fully restored due to the quark masses. The phase transition between phases 
1 and 3 is a crossover, like it is at 9 = 0. The crossover line is defined by the inflection points d 2 (a) /dT 2 = 0. 

At high temperature also the CP-violating phase disappears. This is consistent with the fact that at high 
temperature instanton effects become exponentially suppressed j2Z|- The CP-violating phase is after all realized 
due to the instanton induced interaction. The maximum value of the critical temperature as function of c is 219 
McV. 

We have verified that also for nonzero temperature the minimum of the effective potential is at 9 = 0, which 
means the Vafa-Witten theorem continues to hold in the NJL model at nonzero temperature. The same applies to 
baryon and isospin chemical potential. We also have checked whether there are any local minima in the effective 
potential at nonzero temperature and density, but we found none. 
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Now wc will briefly consider nonzero baryon chemical potential [1b = \i u + [id, where [i u .d denote the u, d quark 
chemical potentials. The ((1b, c) phase diagram is displayed in Fig.[8]for a restricted range of fis values. The same 
phases occur as in the (T, c) phase diagram, but now the phase transition to the (almost) chiral symmetry restored 
phase is of first order, like for 9 = 0. Furthermore, the first-order phase transition has a small c dependence. As 
always, the phase transition from phase 1 to phase 2 is of second order. 



VII. NONZERO ISOSPIN CHEMICAL POTENTIAL 



In quark matter systems equilibrium and neutrality conditions can require that fi u ^ fid- Son and Stephanov 
[38[ observed that charged pion condensation can occur for nonzero isospin chemical potential \±i = \i u — [id- At 
9 = this second order phase transition between the ordinary phase of broken chiral symmetry ((er) ^= 0) to the 
pion condensed phase (which also breaks chiral symmetry) occurs when fii equals the vacuum pion mass. In this 
subsection we address this issue at 9 = n. 

In Fig. [9] we show the phase diagram of the NJL model in the (///, c)-plane, for m u = m,d = 6 McV. The solid 
line indicates a first-order phase transition, the dashed lines indicate second-order phase transitions. The four 
phases are characterized as follows 
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FIG. 9: The (/Ltj, c) phase diagram of the NJL model at 9 = n. 



Phase H] is a novel phase characteristic of 9 = it. This phase also has a small nonzero (cr)-condensate (not 
indicated), due to the explicit breaking by the quark masses. 

For c < c cr it a nonzero (7r ± )-condcnsate exists above a certain fij value. Like at 9 = the second order phase 
transition turns out to be at /i/ = m n , where is the vacuum pion mass. In addition, there is a second phase 
transition, a first-order phase transition, at larger fij, where charged pion condensation makes way for charged ao 
condensation. For c > c cr ;t no nonzero (7r ± )-condensate exists, only nonzero (aj). The phase transition between 
phases [3] and 2] is of first-order. The question arises what determines the value of fii at this phase transition to 
charged meson condensation? To answer this question, the meson masses need to be calculated. 
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A. The c-dependence of the meson masses and mixing 

As said, at 9 = charged pion condensation occurs when /x/ is larger or equal to the vacuum (/zj = 0) pion 
mass. For the NJL model this has been studied extensively in Refs. [35|, bMJ SO, 5jJ ■ This condition is independent 
of c. To see what happens in the 9 = tt case, we calculate the c-dependence of the meson masses, with = 0. 
The results are shown in Fig. [TO] and in Fig. [TT] the c-dependence of the mixing is shown. Clearly, at 9 = tt the 
situation quite different from 9 = 0. 
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FIG. 10: The c-dependence of the meson masses at 9 — tt. The masses are calculated in the RPA. 




FIG. 11: The c-dependence of the mixing-angle of the mesons at 6 — tt. The mixing is calculated in the RPA. 



At c = (no instanton interactions) the rj and 7r masses are equal and also the a and ao masses. This follows 
from the symmetry of the Lagrangian at c = 0, which has a U(2)l<8> U(2)/j-symmetry that is spontaneously broken 
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by the chiral condensate (ignoring the explicit breaking by the quark masses) to V(2)y- This means there are 
four (pseudo-)Goldstone bosons, with the same (small) masses: the 77 and ir mesons. The instanton interactions 
remove the degeneracy for c^O. 

When c > c cr it, i.e. when (77) ^ 0, a complication arises: the mass eigenstates are not CP or P eigenstates 
any longer. The occurrence of the 77 condensate results in mixing of the c-particle with its parity partner, the 77- 
particle. Similarly, the pions mix with their parity partners, the do's. The mixing is to be expected because when 
the ground-state is not CP-conserving, there is no need for the excitations, i.e. the mesons, to be CP eigenstates 
or states of definite parity in case of charged mesons. 

The mass eigenstates, denoted with a tilde, arc defined in the following way 

\a) = cos9 v \a) + sin9 v \rf) , 

\fj) = cos 6^ \r]) - sm9 v \a) , 

\a ) = cos W I ao) + sin 8 n \ ir) , 

\tt) = cos^ |7r) — sin&n- |ao) , (25) 

where 9 V and 9^ are the mixing angles. The states on the r.h.s. are the usual states of definite parity. 

The calculation of the mixing and the resulting masses is similar to the mixing of 770 and 77s in the three-flavor 
NJL-model, which was discussed in great detail in Ref. [42| using the random phase approximation (RPA), the 
approach we will also employ here. As a side remark, we mention that we have also calculated the curvature of 
the effective potential at the minimum which should approximately be proportional to the masses, giving indeed 
very similar results 1 . 

When c < c cr i t no mixing takes place and the tilde fields are equal to their counterparts without tilde. When 
c > c cr it mixing occurs. The mixing between 77 and a increases rapidly as c increases, reaching a maximum at 
c = 0.09, where a is almost completely 77 and vice versa. For larger c the mixing however decreases rapidly again, 
so that when c = \, a (77) is again equal to a (77). 

The mixing between ao and the pions behaves differently, here the mixing angle increases rapidly to become 
90° at c = 5, i.e. 7r becomes ao and vice versa. 

Now we turn to the behavior of the tilde-meson masses, which also display unusual features as function of c. 
When c < c cr ;t, the 7r masses are constant, and equal to the ordinary pion masses. Furthermore, the 77 mass 
decreases with increasing c. This is peculiar to 9 = tt, because at 9 = the 77 mass increases with increasing c. 
The 77 mass has its lowest, nonzero value at c cr i t . This is in contrast to three-flavor lowest order chiral perturbation 
theory [ll|, where the 77 mass (in Ref. actually the tt° mass using the primed theory) vanishes at the phase 
transition. Finally, also the masses of the a 's and tr decrease slightly with increasing c. 

When c > c cr it the c-dependence of the masses changes dramatically. The tt mass now decreases monotonically 
with increasing c, whereas the ao and a mass both increase monotonically to infinity towards c = i. The latter 
can be understood, because when c=\,G\ equals G2 , which as mentioned means for 9 = tt that there arc no 7r 
and a mesons in the spectrum. 

Another striking feature is that the 77 mass rises until it almost reaches the a mass, after which it remains 
approximately constant. The behavior of the a mass is opposite, first it is almost constant and when it becomes 
almost equal to the mass it increased to infinity. The masses of a and 77 cannot cross when there are interactions 
that mix the two states, which is similar to level repulsion in quantum mechanics. The point where both masses 
are almost equal corresponds to a mass that is twice the constituent quark mass. This forms the threshold to 
decay into two quarks, which makes one of the two mesons unstable when c > c cr it- 

Now we turn again to the original question concerning the charged meson condensation phase transition. From 
the calculation of the masses of the tilde-mesons, we infer that the condition for charged meson condensation at 
9 = tt is /J/ > m^(c). For c < c cr it, the phase transition takes place when /i/ equals it does at 9 = 0. 

For c > c cr it it takes place at the mass of tt, which is now a mixed state of 7r and ao. At c = | this means at the 
mass of the ao- The latter observation is in agreement with a result of Ref. [l3[, where the (/// , 9) phase diagram 
of degenerate two-flavor chiral perturbation theory is investigated to lowest order at effectively c — i (due to the 
absence of the 77 meson). There it is observed that charged pion condensation occurs when /jj is equal to the 
^-dependent pion mass m- K {9). In Ref. [l3[ all ^-dependence resides in the mass matrix, with both quarks having 



1 This method was also used in Ref. |39| to calculate the pion mass in order to check the condition fj,f > for charged pion 
condensation. A discrepancy of approximately 25% was obtained. 
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a ^-dependent mass. Hence, their (^-dependent pion field corresponds to what we call 7r', which at 8 = tt is the ao 
field in the original, unprimcd theory, leading to an agreement with our finding. 

The second phase transition for c < c cl -it from the charged pion condensed phase to the charged ao condensed 
phase does not correspond to /jj being equal to a meson mass calculated at /jj = 0. Although we have not found 
a condition in terms of the calculated masses, /jj at this first-order phase transition follows a line that is the 
smooth continuation of the tt mass in the reg ion c > c cr it to infinity at c = 0. A calculation of the meson masses 
at nonzero /ij, such as performed in Ref. [Uj, might resolve this open issue. 




The ^-dependence of the ground-state of the two-flavor NJL model is investigated in the mean-field approxima- 
tion. The main focus is on the case 8 = 7r, when spontaneous CP-violation is possible. The ^-dependence of the 
theory is found to strongly depend on the strength of the 't Hooft determinant interaction. When the strength 
of this interaction, which is governed by the parameter c, is small or zero, no spontaneous CP-violation takes 
place at 8 = tt. The low-energy physics is then almost the same as at 8 = 0. At larger c however, spontaneous 
CP-violation docs take place at 8 = tt. So the phenomenon of spontaneous CP-violation is governed by the 
't Hooft determinant interaction, which describes the effect of instantons in the effective theory. The question 
whether c is sufficiently large for CP-violation to occur at 8 = tt depends on the quark masses. In other words, 
spontaneous CP-violation requires instantons, but its actual realization depends on the size of their contribution 
w.r.t. the quark masses. This is also expected to be the case in QCD, where it can be phrased in terms of the 
low-energy theorem identity ^2 q 2im q ((77.5(7} = —Nj (g 2 FF) /8tt 2 (cf. e.g. Ref. [IH), which relates (77) to the first 
derivative of the effective potential w.r.t. 8. Depending on m q the coupling constant g needs to be sufficiently 
large for spontaneous CP-violation to take place. Or in other words, the energy needs to be sufficiently low. The 
latter observation is in agreement with the disappearance of the CP-violation at temperatures above a certain 
critical temperature or density. Therefore, we conclude that spontaneous CP-violation in the strong interaction 
is an inherently low-energy phenomenon. 

We have checked that the Vafa-Witten theorem holds in the NJL model also at finite temperature and density 
and found that no local minima arise, indicating the absence of meta-stable CP- violating states in the NJL model. 
We have confirmed several previous results that were obtained in two-flavor chiral perturbation theory. We found 
(in accordance with the results of Ref. 0]) that two-flavor lowest-order chiral perturbation theory with a/N — > 00 
is in general not rich enough to yield results that one might expect to hold in QCD too. It leads for instance 
to the conclusion that only for m u = md spontaneous CP violation occurs, without a critical strength of the 
instanton induced interaction. In contrast, the phase diagram of the two-flavor NJL model is very similar to that 
of three-flavor chiral perturbation theory |ll| , where spontaneous CP violation arises for specific ranges of quark 
masses. 

We also found that the presence of a nonzero 77-condensate has a strong effect on the c-dependence of the 
meson masses and gives rise to mixing among the states of definite parity, as expected when CP invariance is 
not a symmetry anymore. As a result, the pions mix with their parity partners, the ao's, and the 77 meson mixes 
with its parity partner, the a meson. Unlike the mixing discussed as a function of 8 which is just a matter of 
consistently naming the states in order to be able to compare to results obtained with negative quark masses 
and which does not affect physical results, the mixing as function of c does change the physics. For instance, 
the condition for charged pion condensation at nonzero isospin chemical potential becomes modified. At 8 = tt 
for c < c cr it, a second-order phase transition takes place when /jj equals m„ just as at 8 = found by Son and 
Stcphanov. However, we find that for c > c cr it it becomes a first-order phase transition to a novel phase of charged 
do condensation that takes place at the mass of tt, which is a mixed state of 7r and a$. At c = h it is entirely ao- 
Charged ao condensation also arises for c < c cr it and fii > m v , but it appears there is no condition in terms of 
vacuum meson masses for this second phase transition. 

We expect the presented two-flavor NJL model results to remain valid in the case of three flavors and when 
going beyond the mean-field approximation, but this remains to be studied. It would be very interesting if the 
results could in the future be compared to lattice QCD results on the low-energy physics at 8 = tt. 




VIII. CONCLUSION AND DISCUSSION 



15 



Acknowledgments 

We would like to thank Harmon Warringa for kindly sharing with us his code to calculate the effective potential. 
We also thank the members of the theory group at the VU, in particular Wilco den Dunnen and Erik Wessels, 
for fruitful discussions. 



[1] V. Baluni, Phys. Rev. D 19, 2227 (1979). 

[2] R. Crewther, P. di Vecchia, G. Veneziano, and E. Witten, Phys. Lett. 88B, 123 (1979). 
[3] K. Kawarabayashi and N. Ohta, Prog. Theor. Phys. 66, 1789 (1981). 
[4] C. Vafa and E. Witten, Phys. Rev. Lett. 53, 535 (1984). 
[5] R. Dashen, Phys. Rev. D 3, 1879 (1971). 
[6] E. Witten, Ann. Phys. 128, 363 (1980). 

[7] P. di Vecchia and G. Veneziano, Nucl. Phys. B 171, 253 (1980). 

[8] A. V. Smilga, Phys. Rev. D 59, 114021 (1999). 

[9] M. H. G. Tytgat, Phys. Rev. D 61, 114009 (2000). 
[10] G. Akemann, J. T. Lenaghan, and K. Splittorff, Phys. Rev. D 65, 085015 (2002). 
[11] M. Creutz, Phys. Rev. Lett. 92, 201601 (2004). 
[12] M. Metlitski and A. Zhitnitsky, Nucl. Phys. B 731, 309 (2005). 
[13] M. Metlitski and A. Zhitnitsky, Phys. Lett. B 633, 721 (2006). 
[14] T. Fujihara, T. Inagaki, and D. Kimura, Prog. Theor. Phys. 117, 139 (2007). 
[15] G. 't Hooft, Phys. Rev. D 14, 3432 (1976). 
[16] G. 't Hooft, Phys. Rep. 142, 357 (1986). 
[17] T. D. Cohen, Phys. Rev. D 64, 047704 (2001). 
[18] M. B. Einhorn and J. Wudka, Phys. Rev. D 67, 045004 (2003). 
[19] T. D. Lee, Phys. Rev. D 8, 1226 (1973). 
[20] P. D. Morley and I. A. Schmidt, Z. Phys. C 26, 627 (1985). 

[21] D. Kharzeev, R. D. Pisarski, and M. H. G. Tytgat, Phys. Rev. Lett. 81, 512 (1998). 

[22] K. Buckley, T. Fugleberg, and A. Zhitnitsky, Phys. Rev. Lett. 84, 4814 (2000). 

[23] D. Kharzeev and R. D. Pisarski, Phys. Rev. D 61, 111901 (2000). 

[24] S. A. Voloshin, Phys. Rev. C 70, 057901 (2004). 

[25] D. Kharzeev and A. Zhitnitsky, Nucl. Phys. A 797, 67 (2007). 

[26] D. E. Kharzeev, L. D. McLerran, and H. J. Warringa, Nucl. Phys. A 803, 227 (2008). 

[27] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 (1977). 

[28] R. D. Peccei and H. R. Quinn, Phys. Rev. D 16, 1791 (1977). 

[29] F. Wilczek, Phys. Rev. Lett. 40, 279 (1978). 

[30] Y. Nambu and G. Jona-Lasinio, Phys Rev. 122, 345 (1961). 

[31] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 124, 246 (1961). 

[32] M. Buballa, Phys. Rept. 407, 205 (2005). 

[33] K. Fujikawa, Phys. Rev. Lett. 42, 1195 (1979). 

[34] M. Frank, M. Buballa, and M. Oertel, Phys. Lett. B 562, 221 (2003). 

[35] H. J. Warringa, D. Boer, and J. O. Andersen, Phys. Rev. D 72, 014015 (2005). 

[36] H. J. Warringa, hep-ph/0606063. 

[37] D. J. Gross, R. D. Pisarski, and L. G. Yaffe, Rev. Mod. Phys. 53, 43 (1981). 
[38] D. T. Son and M. A. Stephanov, Phys. Rev. Lett. 86, 592 (2001). 

[39] A. Barducci, R. Casalbuoni, G. Pettini, and L. Ravagli, Phys. Rev. D 69, 096004 (2004). 
[40] A. Barducci, R. Casalbuoni, G. Pettini, and L. Ravagli, Phys. Rev. D 71, 016011 (2005). 
[41] L. He and P. Zhuang, Phys. Lett. B 615, 93 (2005). 
[42] S. Klevansky, Rev. Mod. Phys. 64, 649 (1992). 



